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ABSTRACT 


, INMESH  is  a  computer  program  for  the  generation  of  two- 
dimensional  boundary-fitted  coordinate  systems.  This  general 
program  creates  a  numerical  transformation  which  maps  an  arbi¬ 
trary  region  in  physical  space  onto  a  regular  region  with  a 
uniform  grid  in  computational  space.  The  computational  region 
may  consist  of  one  or  more  segments  joined  together  so  that 
any  topological  complexity  in  the  physical  region,  such  as 
multiple-connectedness,  is  suitably  resolved.  Numerical 
solutions  of  partial  differential  equations  may  be  obtained 
from  this  computational  field  regardless  of  boundary  shape. 
Emphasis  is  placed  on  ease  of  use  and  interactive  graphics  so 
that  a  mesh  can  be  generated  quickly,  viewed  interactively, 
and  regenerated  immediately  with  slight  alterations  until  a 
suitable  mesh  configuration  has  been  obtained.  A  brief  explana¬ 
tion  of  the  method  is  given  along  with  a  detailed  description 
of  the  input  required  to  produce  boundary-fitted  coordinate 
systems.  Several  sample  meshes  are  shown  for  configurations 
involving  foil-like  bodies.  - 
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ADMINISTRATIVE 


INFORMATION^ 


This  research  was  carried  out  under  the  Mathematical  Sciences  Research  Program, 
Program  Element  61153N,  Task  Area  SR-0140301,  Work  Unit  1840-040,  and  was  supported 
by  the  Naval  Sea  Systems  Command,  Code  SEA  05R24. 


INTRODUCTION 

Research  in  the  area  of  Numerical  Mesh  Generation  has  increased  dramatically 
during  the  last  decade.  This  is  due  in  large  part  to  rapid  advances  both  in  numer¬ 
ical  and  mathematical  techniques  and  in  the  large-scale  computing  capabilities  need¬ 
ed  to  support  these  techniques.  The  value  of  boundary-fitted  coordinate  systems  as 
a  tool  for  the  numerical  solution  of  partial  differential  equations  on  arbitrarily- 
shaped  regions  was  illustrated  in  a  1974  paper  by  Thompson, 1*  in  which  a  surface- 
oriented  grid  was  created  about  an  irregularly-shaped  two-dimensional  "rock."  A 
finite-difference  scheme  was  then  used  to  compute  a  simple  flow  solution  around  this 
body.  The  mesh-generating  procedure  consisted  of  numerically  solving  a  system  of 
elliptic  equations  to  produce  a  numerical  transformation  that  mapped  the  two- 
dimensional  physical  region  onto  a  rectangular  computational  region  which  was 


*A  complete  list  of  references  is  given  on  page  29. 
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over-laid  with  a  uniform  square  grid.  This  mapping  was  defined  so  that  boundaries 
of  the  physical  region  were  transformed  into  boundaries  of  the  computational  region. 
The  inverse  image  of  the  computational  mesh  is  the  boundary- fitted  coordinate  system 
in  which  each  boundary  is  coincident  with  a  constant  coordinate  line.  As  with  the 
solution  of  the  transformation,  the  solution  of  the  flow  problem  can  then  take 
place  in  the  computational  plane  (after  a  transformation  of  the  governing  equations) 
without  regard  to  boundary  shape.  Thompson's  study  served  to  initiate  much  of  the 
subsequent  research  into  numerical  mapping  techniques. 

Although  the  original  mapping  presented  by  Thompson  was  a  relatively  simple 
one,  the  basic  procedure  was  applicable  to  many  physical  problems  involving  field 
solutions.  In  computational  fluid  dynamics,  for  example,  f inite-d if ferencce  schemes 
require  a  grid  in  the  region  of  Interest  so  that  the  flow  characteristics  can  be 
found  at  the  nodal  points.  Difficulties  often  arise  in  such  solution  methods  as  a 
result  of  inaccurate  numerical  representation  of  boundary  conditions.  The  bound¬ 
ary  conditions  are  usually  best  represented  when  the  boundaries  themselves  are  coor¬ 
dinate  lines  so  that  no  interpolation  is  necessary  due  to  irregular  cells  or  compu¬ 
tational  "stars."  Numerically  generated  boundary-f Itted  coordinates  showed  promise 
of  becoming  a  very  valuable  tool,  not  only  in  the  field  of  computational  fluid  dynam¬ 
ics  but  also  In  the  areas  of  heat  transfer,  structural  mechanics,  and  others.  Another 
important  feature  is  that  the  method  is  not  restricted  to  two  dimensions  as  are  some 
other  widely-used  mapping  techniques  such  as  complex  variables. 

At  DTNSRDC,  we  are  interested  in  flow  problems  involving  complex  geometries 
that  often  include  deformable  boundaries  whose  shape  is  not  known  in  advance  but  is 
part  of  the  solution,  as  in  free-surface  problems.  Other  methods  had  been  used  to 
simplify  the  computations  necessary  to  solve  these  challenging  problems,  but  boundary- 
fitted  coordinates  offered  an  efficient  way  of  adapting  the  grid  system  to  geom¬ 
etries  that  changed  with  time.  In  add*^tlon,  all  calculations  are  done  in  the  comput¬ 
ational  space  which  remains  invariant  even  though  the  shape  of  the  physical  region 
may  change.  This  feature  makes  the  method  especially  attractive  for  nonlinear  water 
wave  problems  since  the  solution  of  these  problems  must  be  obtained  at  boundaries 
whose  shape  is  not  known  in  advance.  Finite-difference  schemes  using  boundary- f itted 
coordinates  have  been  used  successfully  for  free-surface  problems  involving  a 
translating  circular  cylinder,  a  moving  transom  stern,  and  breaking  wave  develop¬ 
ment  .  2 .  ^ 
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In  applying  the  boundary-fitted  coordinate/f inite-d if ference  method  to  fluid 
flow  problems,  we  found  that  the  original  type  of  configuration,  which  mapped  the 
entire  physical  region  to  a  single  rectangle,  was  not  suitable  for  all  of  our  geom¬ 
etries.  To  avoid  rewriting  the  mesh  generating  routines  for  each  new  geometry,  we 
needed  a  general  program  that  could  be  used  for  more  complex  situations.  A  decision 
was  made  to  construct  a  grid-generating  program  that  could  easily  be  applied  to  a 
wide  variety  of  problems.  The  program  had  to  be  one  that  was  not  restricted  to  any 
predetermined  mapping  configuration,  but  that  could  be  designed  to  accommodate  even 
the  most  complicated  geometric  domains,  including  those  with  multiple-connectedness. 
This  versatility  can  be  achieved  through  the  use  of  segmented  computational  regions 
tailored  to  handle  the  topological  complexities  of  each  individual  problem.  The 
proper  choice  of  mapping  configuration  is  particularly  important  for  problems  con¬ 
taining  two  or  more  boundaries  of  different  shape,  especially  if  these  boundaries 
vary  with  time.  The  program  also  had  to  be  modular  in  design  so  that  the  subroutines 
which  linked  the  computational  grid  segments  together  could  be  used  not  only  for  the 
mesh  generation,  but  also  for  numerical  solution  of  the  flow  equations.  INMESH  is 
the  result  of  this  effort. 


DESCRIPTION  OF  THE  METHOD 

The  first  computer  program  written  at  DTNSRDC  to  generate  finite-difference 
meshes  for  arbitrary  doubly-connected  two-dimensional  regions  was  called  NUMESH  and 
was  documented  in  a  1977  report.^  NUMESH  could  produce  only  two  different  grid 
c</n f igurat ions .  The  basic  method  of  numerical  mesh  generation  was  extended  to  three- 
dimensional  regions  in  a  1978  report. ^  The  latter  study  produced  some  simple  three- 
dimensional  meshes  about  ship-like  bodies.  The  idea  of  using  a  set  of  rectangular 
grid  segments  to  "build"  a  suitable  computational  region  for  two-dimensional  problems 
was  the  basis  of  a  new  program,  INMESH,  written  in  1982.^  The  current  work  gives  a 
description  of  the  latest  version  of  INMKSH,  along  with  information  needed  for  execu- 


t 


t  the  program. 


Several  new  features  have  been  added  and  the  input  has  been 


The  method  used  to  map  the  physical  region  onto  the  segmented  computational  grid 
is  based  on  a  system  of  either  Laplace  or  Poisson  equations.  The  desired  transforma¬ 
tion  from  (x,y)-space  to  (Cjn)-space  is  represented  by 

C  =  C(x,y)  and  n  =  n(x,y)  (1) 

The  transformations  are  required  to  be  solutions  of  the  equations 

Cxx  +  ^yy 

and  (2) 

hxx  +  ilyy  =  Q(5»n) 


subject  to  appropriate  boundary  conditions.  The  source  functions  P  and  Q  may  be 
zero  (a  Laplace  generating  system)  or  nonzero  (a  Poisson  system)  in  order  to  influ¬ 
ence  the  spacing  of  coordinate  lines  in  the  physical  region.  At  present,  spacing  of 
both  ^  and  n  lines  in  the  physical  region  can  be  made  to  assume,  spacing  similar 
to  that  of  points  on  the  boundaries.  This  feature,  along  with  the  proper  choice  of 
mapping  configuration,  should  provide  enough  coordinate  system  control  for  most 
geometries . 

For  computational  purposes.  Equations  (2)  are  transformed  to  (^,ri)-space  by 
i  nt>'rchanging  dependent  and  independent  variables  to  yield 


whe  tr- 


ax^^  -  2  8x^^  4-  +  j2(Px^  +  Qx^)  =  0 

*  ^yrm  ^  +  Qyn)  =  o 


(3) 


2  2 

a  =  Xt^  >  y^  e  =  x^x^  +  y^y^ 

(4) 

2  2 

y  =  X^  +  y^  J  =  x^y^  -  x^y^ 
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Each  derivative  in  Equations  (3)  and  (4)  is  replaced  by  the  appropriate  central 
difference  approximation,  and  the  resulting  system  is  solved  using  successive  over¬ 
relaxation.  A  more  complete  discussion  of  this  solution  technique  can  be  found  in 
previous  work.^ 

The  computational  region  comprises  rectangular  grid  segments  which  are  mapped 
onto  sections  of  grid  in  the  physical  region.  The  sides  of  the  grid  segments  in 
computational  space  are  joined  together  so  as  to  achieve  the  desired  distribution  of 
corresponding  grid  lines  in  physical  space.  The  sides  of  each  rectangle  are  desig¬ 
nated  as  one  of  two  fundamental  types:  1)  Fixed  (F-type)  and  2)  Matching  (M-type). 

For  the  F-type  side,  the  x-  and  y-coord inates  of  the  grid  points  are  fixed  and  hence 
do  not  change  during  the  computation  of  the  mesh.  An  F-type  side  would  be  found  in 
a  grid  section  adjacent  to  a  boundary  in  physical  space  such  as  a  wall  or  far-field 
bou  dary.  M-type  designations  are  used  to  connect  the  sides  of  different  segments 
or  different  sides  of  the  same  segment.  The  two  sides  that  are  joined  must  both  be 
M-type  and  have  the  same  number  of  points.  An  M-type  connection  has  an  overlap  of 
two  grid  lines.  Unlike  F-type  sides,  the  physical  coordinates  of  M-type  sides  are 
not  determined  in  advance  but  must  be  computed.  This  property  ensures  the  smoothness 
of  coordinate  lines  throughout  the  physical  region.  The  grid  in  the  computational 
region  must  be  defined  so  that  each  side  of  each  segment  is  either  fixed  to  a  phys¬ 
ical  boundary  or  is  joined  to  another  entire  side.  This  requirement  may  mean  that 
many  computational  segments  are  needed  for  a  complex  geometry,  but  the  input  for 
(I’ch  segment  is  simple.  Another  approach  would  have  been  to  allow  the  joining  of 
partial  sides  of  segments.  If  the  computational  segments  were  not  required  to  join 
•'ntire  s.de  to  entire  side,  fewer  of  them  might  be  needed  for  a  particular  geometry, 
hit  the  input  that  designated  the  connections  would  be  more  complicated. 

USE  OK  THE  PROGRAM 

To  create  ,i  coordinate  system  with  INMESH,  the  user  must  first  decide  on  the 
■nipping  configuration  that  best  suits  the  geometry  of  the  problem.  This  step  often 
.iivol.o'S  making  a  sketch  of  tlie  coordinate  system  as  it  will  appear  in  the  physical 
p.',an>’.  Hie  grid  in  the  physical  region  is  then  divided  into  sections  in  s  .ch  a  way 
as  to  ‘:usur>’  t'nat  each  side  of  -each  segment  in  the  compu  ta  t  iona  1  region  is  either 
lixed  to  a  physical  boundary  or  is  joined  to  anoth.-r  entire  side.  Tho  sides  are  de- 
t  in-’d  as  iine  of  the  two  types  (F  or  M),  with  coniu'ctions  spi'ci  fil'd  where  appropriate. 


des  to  be  used  in  the  initialization  of  the  mesh  are  set  along  with  the  coordinate 
ntrol  parameters  for  each  side.  Finally,  the  physical  coordinates  of  the  fixed 
des  are  determined  and  the  program  is  ready  to  be  run. 

All  user-supplied  input  to  the  program  is  given  in  the  form  of  80-column  card 
•iges  with  variables  in  free  format.  The  input  consists  of  values  of  parameters 
ich  1)  define  the  mesh  segments,  2)  designate  the  type  of  each  side,  3)  specify 
des  used  tor  initialization,  4)  establish  the  linkage  between  segment  sides,  and 
provide  coordinate  line  control.  The  input  for  each  segment  side  is  simple--two 
rameters  are  needed  for  connections  and  two  for  control.  The  physical  coordinates 
Che  fixed  sides  must  be  specified  where  necessary.  Also  included  in  the  input 
values  of  parameters  needed  for  the  iterative  solution  of  the  mesh.  The  INMESH 
er's  Guide  in  Appendix  A  gives  a  complete  description  of  the  necessary  input  for 
(MESH,  along  with  a  simple  example. 

In  addition  to  the  printed  output  of  INMESH,  which  echoes  the  input  and  gives  a 
cord  of  the  convergence  error  after  each  iteration,  two  files  can  be  created  for 
iLor  US'-'.  A  plot  file  in  a  format  suitable  for  use  with  the  DTNSRDC  Interactive 
ita  Disfilay  System  (  IDDS)  or  Apollo  Interactive  Data  Display  System  (AIDDS)  plotting 
ickages  can  be  written  as  well  as  a  restart  file  for  use  as  input  to  another  INMESH 
Ml  or  as  the  initial  mesh  for  an  application  program. 

Almost  any  type  of  boundary-fitted  coordinate  system  can  be  created  with  INMESH. 
!'■  (.oordinate  system  that  is  generated  may  be  easily  changed  by  altering  the  number 
■,  r  i  1  segments  and  connections  between  them  in  the  computational  region.  Grid 
les  tr.'  added  to  or  removed  from  specific  areas  by  changing  the  number  of  points 
'  tue  appropriate  segments.  The  coordinate  lines  in  the  physical  plane  are  control- 
‘d  b'/  'varying  the  number  of  grid  lines  in  particular  segments  and  changing  the 
Minition  of  the  functions  P  and  Q  in  the  generating  equations.  The  flexibility 
jvided  by  INMESH  is  especially  useful  for  multi-body  problems  in  which  different 
si'  on f  igura t  ions  are  needed  in  different  parts  of  the  physical  region. 

Die  current  FORTRAN  version  of  TNMESH  is  being  run  on  an  Apollo  DOMAIN  Scien- 
' 1 ,  '.ng i nee r i ng  Workstation  (SEW).  The  SEW  also  supports  interactive  graphics  so 
i  It  a  mesh  can  be  generated  quickly  and  viewed  immediately.  If  changes  are  neces- 
iry,  the  mesh  can  be  regenerated  from  new  initial  data  or  from  the  restart  file 
'oduced  by  the  previous  run  and  modification  data.  This  process  is  continued 
Uil  a  mesh  possessing  the  properties  necessary  for  a  particlar  application  is 
"tiiued.  The  mesh  would  then  be  transferred  to  a  mainframe  computer  for  subsequent 
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TYPICAL  MESH  SEGMENT 


NORTH 


SOUTH 


The  spacing  of  both  C  and  n  lines  within  an  individual  segment  is  influenced 
y  control  functions  defined  on  each  of  the  segment  sides.  These  control  functions 
re  based  on  the  boundary  point  distribution  of  fixed  sides  of  segments  in  the  mesh, 
ny  fixed  side  may  be  used  for  the  control  function  definition  of  that  or  any  other 
egment  side  as  long  as  the  number  of  points  agree.  With  this  type  of  control,  the 
>  lative  spacing  of  the  coordinate  lines  tends  to  be  the  same  as  that  of  the  points 
n  the  controlling  fixed  sides.  The  parameters  KSEG  and  KSID  on  east  and  west  sides 
pecify  the  controlling  segment  sides  for  the  n  lines.  Control  functions  for 
lie  n  lines  in  the  interior  of  a  mesh  segment  are  then  determined  by  linear  interpo- 
atinn  between  the  control  functions  defined  on  the  east  and  west  sides.  The  KSEG 
mi  KSII)  parameters  for  the  north  and  south  sides  control  the  5  lines  in  a  similar 
lai  ner.  If  KSEG  =  KSID  =  0,  the  control  function  is  assumed  to  be  based  on  a  fixed 
.  ide  having  a  uniform  distribution  of  points. 

Type  (4)  input:  Boundary / ini t ia 1 izat ion  data  in  the  form  of  (x,y)  pairs  is  re- 
(uired  for  all  fixed  sides.  If  NTYP  7*  0  in  the  side  connection  data  card,  MSEG  = 

■IS  ID  =  0  and  boundary  data  must  be  supplied  next  in  the  sequence.  NTYP  >  0  indicates 
Chat  the  boundary  data  is  given  in  normal  order,  i.e.  south  to  north  for  the  east  and 
■ifist  sides  and  west  to  east  for  the  north  and  south  sides.  NTYP  <  0  indicates  that 
the  boundary  data  is  given  in  reverse  order  and  is  to  be  converted  to  normal  order 
by  Che  program  before  any  computation  is  done.  Fixed  sides  used  for  initialization 
grid  segments  are  designated  by  NTYP  •  If  two  opposite  sides  of  a  segment 

ist  and  west  or  north  and  south)  are  specified  as  initialization  sides,  linear 
1  n  t-' rpo  1  a  t  ion  is  used  to  distribute  points  throughout  the  interior  of  the  segment. 
Ivich  segment  with  one  or  more  fixed  sides  should  have  an  initialization  side  or 
T  i  :li>3  sp^^c  i  f  ied  . 


;,x  amp  1  ■* 

Figure  A.l  shows  the  segments  of  the  computational  space  for  a  grid,  Mesh  1, 
iboiit  a  simple  body  located  within  a  square  outer  boundary.  Also  shown  in  Figure 
A.l  is  a  typical  mesh  segment  with  sides  labeled  according  to  its  orientation. 

Table  A.l  gives  the  input  necessary  for  the  creation  of  Mesh  1  illustrated  in  Figure 
A. 2.  An  exploded  view  of  the  physical  space  for  this  configuration  is  included  in 
Figuri'  A.l  for  comparison  with  the  corresponding  computational  space  seen  in  Figure 


Type  (3)  input:  The  side  connection/control  function  parameters,  in  order,  are 
follows : 

NSID,  NTYP,  MSEC,  MSID,  KSEG,  KSID 


Parameter  Description 

NSID  Grid  segment  side  identifier  for  current  segment: 

NSID  =  1,  East/Right  (i  =  il) 

=  2,  North/Top  (j  =  jl) 

=  3,  West/Left  (i  =  1) 

=  4,  South/Bottom  (j  =  1) 

NTYP  Side  type  identifier: 

NTYP  =  0,  Matching 

=  1 ,  Fixed  (initialization  side) 

=  ^2,  Fixed 

(if  NTYP  <  0,  input  is  in  reverse  order) 

MSEG  Matching  segment  identifier: 

MSEG  =  0,  No  matching  (NTYP  i  0) 

=  NSEG  of  Matching  segment  (NTYP  =  0) 
MSID  Matching  side  identifier: 

MSID  =  0,  No  matching  (NTYP  f  0) 

=  NSID  of  Matching  side  (NTYP  =  0) 

(if  MSID  =  -NSID,  points  matched  in  reverse  order) 
KSEG  Control  segment  identifier: 

KSEG  =  0,  Control  based  on  equal  spacing 
=  NSEG  of  controlling  segment 
KSID  Control  side  identifier; 

KSID  =  0,  Control  based  on  equal  spacing 
=  NSID  of  controlling  side 
(if  KSID  =  -NSID,  points  used  in  reverse  order) 
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The  remaining  card  images  supply  basic  information  for  the  individual  grid 
segments . 

(2) :  Segment  definition 

(3) :  Side  connections /control  function  definition 

(4) :  Boundary/initialization  data 

(5) :  Control  function  modification 

INITIAL  RUN 

For  an  initial  run  (INT  =  0),  a  sequence  of  card  images,  having  the  form  shown 
in  the  following  table,  must  be  given. 


Description 


SOR  parameters 

Segment  definition  (FIRST  SEGMENT) 

Side  connection/control  function  definition 
Boundary  data  (if  needed) 

Side  connection/control  function  definition 
Boundary  data  (if  needed) 

Side  connection/control  function  definition 
Boundary  data  (if  needed) 

Side  connection/control  function  definition 
Boundary  data  (if  needed) 

Segment  definition  (SECOND  SEGMENT) 

Side  connection/control  function  definition 
Boundary  data  (if  needed) 


\ 

Type  (2)  input:  The  segment  definition  parameters,  in  order,  are  as  follows: 
NSEG,  IL,  JL 


Parameter  Description 

NSEG  Segment  number 

IL  Number  of  ^  (vertical)  lines  in  current  segment 

Number  of  h  (horizontal)  lines  in  current  segment 


JL 


INMESH  USER'S  GUIDE 


Input  needed  for  INMESH  falls  into  five  main  categories;  (1)  parameters  for 
the  successive  overrelaxation  (SOR)  scheme,  (2)  parameters  that  define  the  size  of 
each  segment,  (3)  parameters  that  specify  the  way  in  which  the  grid  lines  are  match¬ 
ed  and  controlled,  (4)  physical  coordinates  of  the  boundaries,  and  (5)  parameters  to 
modify  the  control  functions.  On  initial  runs,  the  first  four  types  of  input  are 
required  with  designated  fixed  sides  providing  an  initial  guess.  On  subsequent  or 
restart  runs,  in  which  a  previously  generated  mesh  is  used  as  the  initial  guess, 
only  the  iteration  scheme  parameters  need  be  given  although  tlie  control  function 
definition  may  be  modified.  All  data  is  provided  in  the  form  of  80-column  card 
images  with  variables  in  free  format. 

Type  (1)  input:  The  first,  card  image  supplies  the  SOR  parameters  for  both  ini¬ 
tial  and  restart  runs.  These  parameters,  in  order,  are  as  follows: 

El,  E2,  MAX,  RF,  INT 

Description 

Minimum  (x|  and  jyj  included  in  convergence  test; 
typically  El  is  about  IZ  of  max(|x|,|y|) 

Maximum  relative  change  allowed  for  convergence: 
typically  0.001  <  E2  <  0.1 

Maximum  number  of  iterations  allowed  in  one  run 

Relaxation  factor: 

typically  1.2  ^  RF  <  1.7 

Initialization  parameter: 

INT  =  0,  Initial  run 

=  1,  Restart  run 


Parameter 

El 

E2 

MAX 

RF 

INT 


X 

Figure  8  -  Free-Surface  Development  for  iranslating  Hydrofoil 


CONCLUSION 

A  computer  program  (INMESH)  for  the  generation  of  boundary-fitted  coordinate 
systems  has  been  written  which  makes  use  of  segmented  computational  regions.  This 
program  is  not  restricted  to  certain  predetermined  mapping  conf igurati^ns  and  there¬ 
fore  can  be  used  to  produce  a  virtually  unlimited  number  of  different  grid  systems. 
The  user  may  design  a  computational  region  suited  to  a  particular  geometry  which 
makes  INMESH  a  useful  tool  in  the  creation  of  curvilinear  coordinate  systems.  The 
current  version  of  INMESH  represents  a  considerable  improvement  on  the  original 
program  with  the  inclusion  of  coordinate  line  control  capability.  INMESH  is  opera¬ 
tional  on  an  Apollo  DOMAIN  Scientific/Engineering  Workstation  that  can  be  used  to 
generate  and  view  meshes  interactively.  The  use  of  INMESH  was  illustrated  by  a 
series  of  sample  grids  and  by  an  application  demonstrating  that  a  segmented  computa¬ 
tional  region  may  be  used  successfully  in  the  numerical  solution  of  a  fluid  dynamics 
p  rob lem. 
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Figure  7b  -  Physical  Region 

Figure  7  -  Example  Grid  for  Foil  in  Rectangular  Region 
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Figure  6  -  Example  3:  Close  View  of  tlie  Confluence  Regions 

on  the  Multibody  Foil 


EXAMPLE  A 

Both  the  inner  and  outer  boundaries  must  often  be  considered  in  choosing  a  grid 
system  for  a  particular  topology.  As  the  final  example,  consider  a  rectangular 
region  containing  a  foil-like  body.  A  useful  grid  for  this  geometry  has  coordinate 
lines  wrapping  around  the  body  yet  has  a  more  rectangular  nature  near  the  outer 
boundary.  The  computational  and  physical  regions  for  this  mapping  configuration  are 
shcawn  in  Figure  7.  In  Figure  7(b),  note  the  irn^gular  five-sided  cells  that  appear 
both  above  and  below  the  foil  where  the  grid  system  changes  character. 

A  coordinate  system  similar  to  the  one  illustrated  in  Figure  7  has  been  used  as 
the  initial  mesh  for  a  preliminary  calculation  of  the  potential  flow  resulting  from 
a  hydrofoil  moving  beneath  a  f ree-sur face .  A  finite-difference  technique  was  used 
along  with  a  marcViing  scheme  for  the  ailvancemenl  of  the  t  ime-d  ependiMit  f  e'e-surface, 
r.cjuiring  that  a  new  mesh  be  gcMierated  at  each  Lime  st.p  as  tne  f  re  e -s  u  r  t  ac  e  evolves. 
Details  of  this  method  can  be  found  in  an  earli-M  papor  iiy  Haussling  and  Coleman.^ 
Figure  8  shows  the  development  of  the  free  surface  loi  various  tiinos  t  as  tlie  foil 
is  abruptly  put  in  motion.  This  computation  was  .  isily  pertonned  on  the  segmented 
grid  since  the  same  subroutines  used  to  gen'rato  tlio  mesh  were  also  usmi  for  the 
flow  calculations. 


Figure  4  -  Example  2:  Close  View  of  Foil  with  Slot 


XAMPLE  3 

Multibody  problems  can  also  be  handled  by  INMESH  even  though  the  mapping  con¬ 
figuration  becomes  somewhat  more  complicated  as  illustrated  by  Figure  5.  The 
increase  in  complexity  of  the  computational  region  arises  from  the  increase  in  the 
number  of  discontinuous  coordinate  lines  in  the  physical  region.  This  geometry  is 
modeled  after  a  tri-element  Augmentor-Wing  airfoil  being  studied  for  use  in  short- 
takeof  f-and -land  ing  aircraft. 8  a  close'-up  view  of  the  confluence  region  is  given 
in  Figure  6. 


1  I 


DPAC  IS  given  in  Appendix  B.  To  illustrate  the  influence  of  boundary-spacing  control 
upon  the  mesh,  another  grid  was  generated  for  the  same  region  with  no  coordinate 
line  control  (P  =  Q  =  0)*  This  mesh  is  shown  in  Figure  2.  More  information  about 
control  of  the  coordinate  lines  can  be  found  in  Appendix  A. 
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Figure  2  -  Example  I :  Grid  For  Region  Shown  in  Figure  I  Generated 
without  Coordinate  Line  Control 


EXAMPLE  2 

The  second  example  illustrates  a  C-type  grid  generated  about  a  foil  that  has  a 
notch  or  slot  in  the  upper  surface.  Figure  3  shows  the  computational  and  physical 
regions  for  this  geometry.  As  in  Example  1,  the  boundary-s pac ing  control  parameters 
were  used  to  attract  coordinate  lines  into  the  area  near  the  body.  This  is  a  model 
of  the  cross  section  of  an  aircraft  wing  that  has  slot  on  the  upper  surface  from 
which  air  is  exhausted  to  improve  the  lift  characteristics.  Figure  4  is  a  close-up 
view  of  the  slot. 


use,  or  perhaps  another  mesh,  larger  but  with  similar  propperties,  would  be  produced 
before  transfer.  In  this  way,  INMESH  can  aid  the  researcher  in  obtaining  a  "best" 
mesh  for  a  particular  problem  by  enabling  him  or  her  to  efficiently  produce  a  wide 
variety  of  grids  for  the  same  geometry. 


EXAMPLES  AND  APPLICATIONS 

EXAMPLE  1 

As  a  first  example  of  the  meshes  that  can  be  created  with  INMESH,  consider 
the  4-segment  computational  region  of  Figure  la  and  the  associated  physical  region 
of  Figure  lb.  This  figure  illustrates  a  "C-type"  coordinate  system  about  a  foil¬ 
like  body.  In  this  and  succeeding  examples,  the  segmented  computational  region  will 
be  illustrated  schematically  with  connections  among  the  segments  shown  as  curved  or 
straight  line  segments.  Sides  with  no  connections  given  are  fixed.  These  sketches 
of  the  computational  regions  show  how  linkage  of  grid  segment  sides  determines  the 
nature  of  the  coordinate  system  as  it  appears  in  the  physical  region.  The  number  of 
grid  lines  in  each  segment  must  be  such  that  every  pair  of  matching  sides  has  an 
equal  number  of  grid  points. 

Note  that  four  computational  segments  are  required  for  this  example  because  the 
coordinate  line  that  lies  along  the  body  also  lies  within  the  physical  region. 

Segment  4  is  used  to  link  Segment  1  with  Segment  3  to  form  a  "cut"  in  the  physical 
region.  Although  an  additional  mesh  segment  is  required  if  a  coordinate  line 
"splits"  in  physical  space,  this  approach  allows  all  sides  of  all  segments  to  be  one 
of  two  basic  types  rather  than  three  as  in  an  earlier  version  of  INMESH.  This  is 
an  example  of  a  simple  mesh  that  has  many  applications. 

Coordinate  line  control  was  used  in  this  example  to  attract  grid  points  into 
the  areas  around  and  downstream  from  the  body.  The  spacing  of  points  on  the  bound¬ 
aries  was  adjusted  and  then  the  control  parameters  were  used  to  maintain  this  spacing 
throughout  the  mesh.  This  particular  distribution  of  boundary  points  was  obtained 
by  using  a  curve-fitting  procedure  called  Distribute  Points  on  A  Curve  (DPAC).  DPAC 
runs  on  the  Apollo  SEW  and  uses  B-splines  to  distribute  any  number  of  points  on  an 
arbitrary  curve.  The  distribution  can  be  specified  explicitly  or  it  can  be  gener¬ 
ated  by  DPAC  based  on  user-supplied  input  concerning  the  clustering.  More  detail  on 
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TABLE  A.l  -  SAMPLE  INPUT  FOR  INITIAL  RUN 


Data 

.001  .01 

100  1 

.5  0 

1  8 

11 

1  0 

2  3  0 

0 

2  0 

4  4  0 

0 

3  2 

0  0  0 

0 

-1.0 

-1.40 

-1.0 

-1.17 

-1.0 

-0.96 

-1 .0 

-0.77 

-1.0 

-0.60 

-1.0 

-0.45 

-1.0 

-0.32 

-1.0 

-0,21 

-1.0 

-0.12 

-1.0 

-0,05 

-1.0 

0.0 

4  1 

0  0  0 

0 

-1.0 

-1.4 

-0.9 

-1.4 

-0.8 

-1.4 

-0,7 

-1.4 

-0.6 

-1.4 

-0.5 

-1  .4 

-0.4 

-1 .4 

-0.3 

-1 .4 

2  9 

11 

1  0 

3  3  0 

0 

2  2 

0  0  0 

0 

-0.3 

0.0 

-0.3 

0.0 

-0.2 

-0.1 

-0.1 

-0.2 

1 

0.0 

-0.2 

0.1 

-0.2 

0.2 

-0,1 

0.3 

0.0 

0.3 

0.0 

3  0 

1  1  0 

0 

4  1 

0  0  0 

0 

-0.4 

-1.4 

-0.3 

-1 .4 

-0.2 

-1.4 

-0.1 

-1 .4 

0.0 

-1 .4 

0.1 

-1 .4 

Data  con't 
from  Col.  2 


0.2  -1.4 
0.3  -1.4 
0.4  -1.4 
3  8  11 
1  2  0  0  0  0 
1.0  -1.40 
1.0  -1.17 
l.O  -0.96 
1.0  -0.77 
1 .0  -0.60 
1.0  -0.45 
l.O  -0.32 
1.0  -0.21 
1  .0  -0.12 
1.0  -0.05 
1.0  1.0 

2  0  5  4  0  0 

3  0  2  1  0  0 

4  1  0  0  0  0 
0.3  -1.4 
0.4  -1.4 
0.5  -1.4 
0.6  -1.4 
0.7  -1  .4 
0.8  -1.4 
0.9  -1.4 
1.0  -1.4 

4  8  3 

1  I  0  0  0  0 
-0.3  0.0 
-0.3  0.0 
-0.3  0.0 

2  0  6  4  0  0 


Data  con't 
from  Col.  3 


3  1  0  0  0  0 

0.3  0.0 

0.3  0.0 

0.3  0.0 

4  0  3  2  0  0 

6  8  11 

1  0  7  3  0  0 

2  1  0  0  0  0 
-1.0  1.4 
-0.9  1.4 
-0.8  1.4 
-0.7  1.4 
-0.6  1.4 
-0.5  1.4 
-0.4  1  .4 
-0.3  1.4 

3  2  0  0  0  0 

-1.0  0.0 

-1  .0  0.05 
-l.O  0.12 
-1.0  0.21 
-1.0  0.32 
-1  .0  0.45 
-1.0  0.60 
-1 .0  0.77 
-1.0  0.96 
-1.0  1.17 
-1.0  1.40 

4  0  4  0  0  0 

7  9  11 

1  0  8  3  0  0 

2  1  0  0  0  0 


3 

2 

0  0  0 

0 

-0. 

.4 

1 

.4 

-1 , 

.0 

-0.1 

-0. 

.3 

1 

.4 

-1 

.0 

0.0 

-0. 

,2 

1 

.4 

-1 

.0 

0.1 

-0. 

,  1 

1 

.4 

4 

0 

1  2  0 

0 

0. 

.0 

1 

.4 

5 

8 

3 

0, 

.1 

1 

.4 

1 

2 

0  0  0 

0 

0. 

,2 

1 

.4 

1 

.0 

-0.1 

0. 

.3 

1 

.4 

1 

.0 

0.0 

0. 

.4 

1 

.4 

1 

.0 

0.1 

3 

0 

6 

1 

2 

0 

8  4  0 

0 

4 

2 

0 

0 

Data  con ' t 
from  Col.  4 


-0.3  0.0 
-0.3  0.0 
-0.2  0.1 
-0.1  0.2 
0.0  0.2 
0.1  0.2 
0.2  0.1 
0.3  0.0 
0.3  0.0 
8  8  11 
1  2  0  0  0  0 
1  .0  0.0 
1.0  0.05 
1.0  0.12 
1.0  0.21 
1.0  0.32 
1.0  0.45 
1  .0  0.60 
1.0  0.77 
1.0  0.96 
1.0  1.17 

1  .0  1 .40 

2  1  0  0  0  0 
0.3  1.4 
0.4  1 .4 
0.5  1.4 
0.6  1.4 
0.7  1.4 
0.8  1.4 
0.9  1.4 

1.0  1.4 

3  0  7  1  0  0 

4  0  5  2  0  0 


RESTART  RUN 


For  a  restart  run  (INT  =  1),  only  the  SOR  parameter  card  (1)  is  required.  The 
user  also  has  the  option  of  altering  the  control  functions  of  the  previously  gener¬ 
ated  mesh.  This  is  done  by  supplying  one  control  function  modification  card  for 
each  segment  side  on  which  the  function  is  being  redefined.  Thus  the  complete 
sequence  of  card  images  needed  for  a  restart  run  is  as  follows: 

(1)  SOR  parameters 

(5)  Control  function  modification  (if  needed) 

(5)  Control  function  modification  (if  needed) 

etc. 

Type  (5)  input:  The  control  function  modification  parameters,  in  order,  are  as 
follows : 

NSEG,  NSID,  KSEG,  KSID 

Parameter  Description 

NSEG  Segment  for  which  control  is  being  modified 

NSID  Side  on  which  control  is  being  modified 

KSEG  Controlling  segment  identifier 

KSID  Controlling  side  identifier 

Example 

Input  data  used  for  a  restart  run  with  control  function  modification  is  shown 
in  Table  A. 2.  The  output  for  Mesh  1  was  used  as  a  basis  for  thjis  restart,  which 
produced  Mesh  2  (as  shown  in  Figure  A. 3).  A  comparison  of  Mesh  1  and  Mesh  2  shows 
how  the  control  parameters  can  affect  the  curvature  of  the  coordinate  lines  in  the 
physical  space. 


APPENDIX  B 


DPAC:  DISTRIBUTE  POINTS  ON  A  CURVE 

INMESH  uses  an  elliptic  system  of  partial  differential  equations  (a  Laplace  or 
Poisson  system)  to  generate  numerical  grids.  The  curvilinear  coordinates,  5  and 
n,  are  specified  on  the  boundaries  from  which  the  Cartesian  coordinates,  x  and 
y,  are  generated  in  the  bounded  region.  The  boundary  value  problem  is  then  solved 
on  the  physical  field  with  boundary  conditions  specified  on  curved  boundaries. 

Correct  specification  of  these  boundary  conditions  is  thus  crucial  to  obtaining  the 
desired  mesh.  DPAC  is  a  set  of  numerical  procedures  designed  to  aid  the  user  in 
producing  a  distribution  of  boundary  points  so  that  a  suitable  coordinate  system 
results.  While  the  interactive  curve-fitting  capability  provided  by  DPAC  is  not 
necessary  for  grid  generation  with  INMESH,  it  is  highly  desirable.  The  ability  to 
change  the  boundary  point  distribution  quickly  and  easily  becomes  even  more  important 
as  the  complexity  of  the  coordinate  system  increases. 

DPAC  was  originally  written  for  use  with  the  IDDS  plotting  package  on  the  CDC 
176  at  DTNSRDC.  The  procedures  have  been  converted  to  run  on  the  Apollo  DOMAIN 
Sc ientif ic /Engineering  Workstation  in  conjunction  with  the  AIDDS  plotting  system. 

To  the  user,  DPAC  appears  as  a  series  of  interactive  commands.  After  the  user  enters 
a  command  with  the  appropriate  parameters,  the  corresponding  procedure  is  executed 
and  the  results  are  displayed  graphically.  The  results,  which  have  been  written  to 
a  specified  file,  can  then  be  accepted  or  else  additional  commands  can  be  entered 
until  satisfactory  results  are  obtained. 

A  report  currently  in  preparation  will  describe  the  use  of  DPAC  in  detail.  The 
following  is  a  brief  summary  of  the  basic  DPAC  commands  used  to  produce  the  boundary 
point  distribution  for  the  examples  in  this  report.* 

DE  F  [ nf i le ]  [  nb  J 

The  DEF  command  produces  a  B-spline  fit  using  "nb"  basis  functions  to  the  curve 
described  in  discrete  form  on  file  "nfile."  The  nb  parameter  determines  the  "stiff¬ 
ness"  of  the  fit.  The  original  points  are  displayed  along  with  the  B-spline  curve 
that  has  been  fit  to  them. 


*For  more  information  on  DPAC,  contact  Mr.  M.  L.  Brabanski,  DTNSRDC,  Code  1843. 
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EVEN  [npl 

The  EVEN  command  distributes  "np”  equally  spaced  points  on  the  curve 
defined  by  DEF.  Both  the  points  and  the  curve  are  displayed.  Also  shown  in 
the  display  is  the  interval  [0,1]  with  an  equivalent  distribution  of  points. 

EVALF  [nfile] 

The  EVALF  command  distributes  points  on  a  curve  defined  by  DEF  in  the 
same  manner  as  a  distribution  on  [0,1]  specified  in  "nfile."  The  display 
shows  the  distribution  on  both  the  curve  and  [0,1]. 

ZONES  [nfile]  [np] 

The  ZONES  command  aids  the  user  in  defining  a  distribution  if  the  desired 
distribution  is  not  known  exactly.  A  distribution  of  "np"  points  on  [0,1]  is 
created  by  specifying  in  "nfile"  subintervals  of  [0,1]  where  the  spacing 
between  points  is  expanding,  contracting,  or  remaining  constant.  The  distri¬ 
bution  thus  defined  is  then  automatically  transferred  to  the  curve  by  EVALF. 
Both  distributions,  on  the  curve  and  on  [0,1],  are  displayed. 
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